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We investigate the condensation phase transitions of conserved-mass aggregation (CA) model on 
weighted scale-free networks (WSFNs). In WSFNs, the weight Wij is assigned to the link between 



the nodes i and j. We consider the symmetric weight given as Wi 



In CA model, the 



mass vm on the randomly chosen node i diffuses to a linked neighbor of i,j, with the rate Tji or an 
unit mass chips off from the node i to j with the rate ioTji. The hopping probability Tji is given as 
Tji = Wjij X^<;> w Hy where the sum runs over the linked neighbors of the node i. On the WSFNs, 
we numerically show that a certain critical a c exists below which CA model undergoes the same 
type of the condensation transitions as those of CA model on regular lattices. However for a > a c , 
the condensation always occurs for any density p and u>. We analytically find a c — (7 — 3)/2 on 
the WSFN with the degree exponent 7. To obtain a c , we analytically derive the scaling behavior 
of the stationary distribution P^° of finding a walker at nodes with degree k, and the probability 
D(k) of finding two walkers simultaneously at the same node with degree k. We find ~ fc a+1_7 
and D(k) ~ fc 2 ' a+1 - 1 ' respectively. With P£° , we also show analytically and numerically that the 
average mass m(k) on a node with degree k scales as k a+1 without any jumps at the maximal degree 
of the network for any p as in the SFNs with 01 = 0. 

PACS numbers: 05.70.Fh,05.40.-a,89.75.Da,89.75.Hc 



I. INTRODUCTION 

A wide variety of mass transport systems ranging from 
traffic flow to polymer gels [EH, B 0, H> B 0, H exhibit 
nonequilibrium condensation phenomena. These include 
basic microscopic dynamics ubiquitous in nature such as 
aggregation, fragmentation and diffusion. The nonequi- 
librium steady states of these systems are classified into 
two types of phases, so-called fluid phase and condensed 
phase. A finite fraction of total particles condenses on a 
single site in the condensed phase. In the fluid phase, the 
particle number on each site fluctuates around the den- 
sity of total particles (p) without the condensation. As 
the rates of these processes vary, the condensation phase 
transitions between the two phases may take place at a 
certain critical density p c . 

One of the simplest mass transport models exhibiting 
the condensation transitions is a conserved-mass aggrega- 
tion (CA) model @, EE El El El ■ CA model evolves via 
diffusion, chipping and aggregation upon contact which 
arise in a variety of phenomena such aspolymer gels , 
the formation of colloidal suspensions [5j , river networks 
[H, 0] and clouds @. In one-dimensional CA model, the 
mass m,i of site i moves either to site i — 1 or to site i + 1 
with unit rate, and then m; — -> and rrii±i — > mm +m,. 
With rate tu, unit mass chips off from site i and moves 
to one of the nearest neighboring sites; — ► m, — 1 
and rrii±i — > vrii±\ + 1. The generalization to higher 
dimensions is straightforward. As total masses are con- 
served, the conserved density p and the rate to determine 
the phase of CA model. The ui = 00 case corresponds to 
the well-known zero r ang e process (ZRP) with a constant 
hopping rate 0, El, Elf- 

The existence of the condensation transitions in CA 



model depends on the symmetry of movement, the con- 
straints of diffusion rate and the underlying network 
structure [n| EE El El- In the symmetric CA (SCA) 
model 0, E3] in which diffusion and chipping direction 
are unbiased, the condensation transitions take place at 
a certain p c . The steady state properties of SCA model 
is exactly described by the mean field theory [n|. The 
single site mass distribution P(m) was shown to undergo 
phase transitions on regular lattices 0. For a fixed uj, as 
p is varied across the critical density p c (j-^>), the behavior 
of P(m) for large m was found to be @ 



P(m) 



e~ m/m ' P<Pc(u), 
mT T p = p c (u), (1) 

m~ T + infinite aggregate p > p c {w). 



Mean field theory predicts p c (u) = y/uT+T — 1 and t = 
5/2 . 

Recently, CA model on unweighted scale-free networks 
(SFNs) with degree distribution P(k) ~ k~ 7 was studied 
to investigate the effect of underlying network structure 
on the condensation transitions [13]. We call networks 
with equal weight on all links unweighted networks. On 
unweighted SFNs, the same type of the condensation 
transitions as those of SCA in regular lattice take place 
for 7 > 3. However for 7 < 3, the condensation always 
occurs for any density p(> 0). It was shown that the exis- 
tence of the transitions is directly related to the diffusive 
capture process on unweighted SFNs [H, HI . 

On the other hand, most real-world networks exhibit 
not only a heterogeneous distribution of degree, but 
also heterogeneous distribution of weig hts Ha El El- 
Weights assigned on links characterize the interaction 



2 



strengths between nodes. There have been various at- 
tempts to understand the underlying mechanism and 
scale- free behaviors of empirically observed weighted net- 
works [2(| • Also there have been attempts to understand 
the effect of heterogeneous weights on various dynam- 
ics such as synchronization, dynamics of random walks, 
transport and percolation, and condensation of zero- 
range process [2l|, [22|, [H[. These studies showed that 
dynamical properties are modified and exhibit non-trivial 
dependence on the strength of weight. In this paper, as 
the generalization of our study on CA model on complex 
networks, we investigate the effect of both heterogeneous 
degree and weight on the condensation phenomena of CA 
model on weighted networks. 

The weight Wji represents the weight to a link from 
the node i to j. In general, the strength Si of the node 
i scales with the degree ki as Sj ~ fcf . The exponent a 
varies with network structures \19L l2(j | . Thus it is natural 
to take the weight Wji as Wji ~ SiSj ~ (kikj) a . 

In this paper, we study the condensation transitions of 
CA on the WSFNs with degree distribution P(k) ~ k" 1 
and the symmetric weight wji = (kikj) a . As in one 
dimension, the diffusion of the whole masses and the 
fragmentation of unit mass occur with the unit rate 
and the rate u>, respectively. In addition, masses move 
from a node i to j with hopping rate proportional to 
Wji j w ji ■ We found that a certain critical a c ex- 
ists below which the condensation transitions take place. 
However for a > a c , the condensation always occurs for 
any density p > 0. To find a c as a function of the degree 
exponent 7, one needs the steady state distribution P£° 
of finding a walker at nodes with degree k on the WSFNs. 
-Pfc° gives the capture probability D(k) with which two 
walkers meet at a node with degree k. We analytically 
derived and D(k), and finally obtained a c = (7— 3)/2. 

This paper is organized as follows. In Sec. II, we dis- 
cuss the condensation transitions of CA model on the 
WSFNs. To verify the existence a c , we investigate the 
steady state property of a single walker and the diffusive 
capture process on the WSFNs in Sec. Ill and IV. We 
discuss the behavior of an average mass m(k) of a node 
with degree k in Sec. V and finally summarize our results 
in Sec. VI. 



II. CA MODEL ON WSFNS WITH 
SYMMETRIC WEIGHTS 

We consider CA model on WSFNs with the weight 
Wij from node j to i defined as iD tJ = (kikj) a . For the 
construction of WSFN, we first construct an unweighted 
static SFN with N nodes and K links [26| . The degree ki 
of a node i is defined as the number of its links connected 
to other nodes. The average degree of a node < k > is 
given as < k >= 2K/N. The degree distribution P{k) 
of SFNs is a power-law distribution P(k) ~ fc~ 7 . In the 
static model [2g], it is desired to use large < k > to con- 
struct fully connected networks. In simulations, we use 



< k >= 4. After then, we assign a weight Wij = (kikj) a 
to the link between node i and j. Thus the hopping prob- 
ability of masses from node i to an i's linked neighbor j 

is T — h a h a I lc a h a — h a I ' V h a V 

10 ± Jl — fitj /Vj / 2^< m > ^m^i ~ "y I £~i<m> Nn- Z^<m> 

denotes the sum over the linked neighbors of node i. 

Each node has an integer number of particles, and the 
mass on a node is defined as the number of particles on 
the node. Initially M particles are randomly distributed 
on N nodes with given conserved density p = M/N. Next 
a node i is chosen at random and one of the following 
events occurs: 

(i) Diffusion : With the unit rate, the whole mass mi 
of node i moves to one of the linked neighbors j with 
probability Z^. Then the aggregation takes place; m,; — ► 
and nij — > rrij + m., . 

(ii) Chipping : With the rate li, unit mass moves to a 
linked neighbor j with the probability Tji, and then the 
aggregation takes place, i.e. m» — ► m$ — 1, rrij — * rrij + 1. 

The u) = 00 case corresponds to ZRP with constant 
chipping rate on WSFNs [I6[ . 

We perform Monte Carlo simulations with random ini- 
tial mass distribution on the WSFNs with 7 = 2.7 and 
4.0. We set w = 1 and the network size N = 10 5 with 

< k >= 4. We measure the single node mass distribution 
P(m) in the steady states. 

In Fig. 1, we plot P(m) for 7 = 4.0 with two differ- 
ent a, a = 0.05 and 1.0. P(m) exhibits quite different 
behavior according to the value of a. For a = 0.05 (Fig. 
1(a)), P{m) decays exponentially without aggregates for 
sufficiently low density p = 0.15. On the other hand, for 
sufficiently high density, p = 3.0, an aggregate forms with 
the power-law decaying background mass distribution. It 
means that the condensation transition takes place at a 
certain critical density p c (> 0). Hence P(m) follows Eq. 
(fl]). Since in unweighted SFNs, i.c a = 0, the conden- 
sation phase transitions take place for 7 > 3 [13j, one 
may expect the condensation transitions for very small 
a. Based on the following steps, we estimate p c and the 
exponent r. 

In the condensed phase, the total density p is written 
as p = p c + poo, where p^ is the density of an aggre- 
gate. Since p is given as p = J ± mP(m)dm, one can 
estimate p c from p c = J" 1 " mP(m)dm, where the upper 
bound m is the cut-off mass at which the background 
distribution terminates. Using this method, we estimate 
p c = 0.218. We estimate the exponent r from the scaling 
plot m T P(m) using P(m) of p = 3.0 (Inset of Fig. 1(a)). 
Since the background distribution docs not change for 
P > Pa we use P(m) of m < m for the scaling plot. We 
estimate t = 2.38(5). 

On the other hand, for a = 1.0 (Fig. 1(b)), P{m) shows 
the complete different behavior. The condensation takes 
place with an exponentially decaying background distri- 
bution for both sufficiently low and high density, p = 0.1 
and 3.0. Therefore we conclude that the condensation 
always occurs for any nonzero density so a system is 
always in the condensed phase without any transitions 
for a = 1.0. The two different behaviors of P{m) for 
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FIG. 1: (Color online) The plot of P(m) for 7 = 4.0 with 
a — 0.05 (a) and a = 1.0 (b). The inset of (a) shows the 
scaling plot m T P(m) with r = 2.38 when p = 3.0. 



a = 0.05 and 1.0 indicate that a crossover a (a c ) should 
exist in the range 0.05 < a < 1.0 for 7 = 4.0. A system 
undergoes the condensation transition for a < a c , while 
the condensation always occurs without the transition for 
a > a c . 

Similarly, for 7 = 2.7, P(m) exhibits the same different 
behavior according to the value of a. The difference from 
the 7 = 4.0 case is that the condensation transitions are 
observed for a negative a. We observe the condensation 
transitions for a = —1.0 (Fig. 2(a)). With the same 
method used in the 7 = 4.0 case, we estimate p c = 0.4 
and t = 2.46(5) respectively. However, for a = —0.05, 
the condensation is observed even for very low density 
p = 0.1, which means that a system is always in the con- 
densed phase for a = —0.05 (Fig. 2(b)). Therefore, the 
crossover a c also exists for 7 = 2.7, but its value is nega- 
tive unlike the 7 = 4.0 case. Together with the results of 
7 = 4.0, we conclude that the crossover a c exists for any 
7(> 2) and a c varies with 7. In what follows, we discuss 
the existence of a c and next the condensation nature for 
a < a c . 

First, the condensation phenomena of CA model on 
WSFNs is similar to that on unweighted SFNs. On un- 
weighted SFNs, the condensation transitions exist for 
7 > 3, while the condensation always occurs for 7 < 3 
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FIG. 2: (Color online) The plot of P(m) for 7 = 2.7 with 
a = —1.0 (a) and a = —0.05 (b). The inset of (a) shows the 
scaling plot m T P(m) with r = 2.46 when p = 3.0. 



[l3j . Hence the crossover 7 is 7 C = 3. Intriguingly, it was 
shown that the existence of transitions is determined by 
the survival probability of a diffusing prey chased by a 
diffusing predator, so-called the lamb-lion problem [271 ]. 
The reason is as what follows. 

In the limit p — > 0, let us assume only an infinite ag- 
gregate exists. With the rate to, the unit mass is chipped 
off and moves around with the unit rate. If the chipped 
mass meets again with the infinite aggregate within a fi- 
nite time interval, then the infinite aggregate is stable 
against the chipping process. On the other hand, if the 
chipped mass and the infinite aggregate does not meet 
again within a finite time interval, then the infinite ag- 
gregate would disappear by repeated chipping processes. 
Therefore, the stability of the infinite aggregate is phys- 
ically related to the capture process in which a diffusing 
lion (infinite aggregate) chases a diffusing lamb (chipped 
mass). For the unweighted SFNs with 7 < 3, it was 
shown that the survival probability S(t) of a lamb decays 
exponentially with finite life time (T)oo [H, HH- How- 
ever, for 7 > 3, S(t) is finite in the thermodynamic limit. 
The behavior of S(t) implies that the condensation tran- 
sition exist for 7 > 3 due to the stable fluid phase in the 
limit p —> 0, but only the condensation exist for 7 < 3. 
As a result, the asymptotic behavior of the survival prob- 
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ability of a lamb in the lamb-lion capture process deter- 
mines the existence of the condensation transitions on 
unweighted SFNs. 

Similarly, on the WSFNs, the existence of the con- 
densation transitions is also expected to depend on the 
survival probability S(t) of a lamb. To see this, let us 
consider two limits, a — > +00 and —00 for a given 7. 
In the limit a — ► +00 , a walker always moves to a node 
with the larger degree. Once a walker reach the hub node 
with the maximal degree, the walker is trapped at the hub 
node forever. As a result, a lion always captures a lamb 
at the hub node within a finite time interval. Hence, S(t) 
should decay exponentially with a finite life time. On the 
other hand, in the limit a — > —00, a walker is forced to 
reach nodes with the minimal degree. Due to the inhomo- 
geneous structure, the nodes with the minial degree are 
connected by nodes with larger degree. Hence, a walker 
cannot escape from one of the nodes with the minimum 
degree in this limit. It means that a lion cannot always 
capture a lamb at some other node, so that S(t) is finite. 

From the behavior of S(t) in the two opposite lim- 
its, there should be a crossover a c . S(t) is finite for 
a < a c and decays to zero for a > a c . For the condensa- 
tion phenomena, one expects no condensation transitions 
(p c = 0) for a > a c due to finite life time of a lamb. 
Instead, the condensation always occurs. On the other 
hand, the condensation transitions occur for a < a c . We 
analytically find a c — (7 — 3)/2 for a given 7 in Sec. IV. 
From a c — (7 — 3)/2, one reads a c = —0.15 for 7 = 2.7 
and a c = 0.5 for 7 = 4 respectively. Our simulation re- 
sults for 7 — 4 and 2.7 confirm the existence of a c and 
also the sign of a c for each 7. 

Next, we discuss the critical behavior of CA model on 
WSFNs. The CA model on any dimensional regular lat- 
tice and unweighted SFNs with 7 > 3 is well described by 
mean- field theory [lj], [l]|. On WSFNs, interestingly, the 
transitions take place even for 7 < 3, which means that 
the transition nature is not affected by the inhomogenc- 
ity of network structure. Since a c diverges for 7 — > 00, 
the critical behavior of CA model on SFNs with a < a c 
should be the same as that on random networks where 
a, i.e. weight, does not have no special meaning due to 
the uniform degree distribution. As a result, one expects 
the mean-field critical behavior of SCA model in regu- 
lar lattice. Our numerical estimates of r, r = 2.38(5) 
for 7 = 2.7 and r = 2.46(5) for 4.0, well agree with the 
mean-field value r = 5/2. Therefore, we conclude that 
the critical behavior of CA model for a < a c on the 
WSFNs belongs to the universality class of SCA model 
in regular lattice. 

In summary, for a fixed 7, there is a crossover weight 
exponent a c . CA model undergoes the same type of con- 
densation transitions as those of SCA model in regular 
lattice for a < a c , while the condensation always takes 
place for nonzero density for a > a c . To find a c as a 
function of the degree exponent 7, one needs the steady 
state distribution P^° of finding a walker at nodes with 
degree k on the WSFNs. In the next section, we derive 



P£° on the WSFNs. In Sec. IV, we study lamb-lion 
capture process on the WSFNs and finally find a c using 

poo 
r k ■ 

III. WALKS ON WSFNS WITH SYMMETRIC 
WEIGHTS 

We consider a single walker on weighted networks with 
the weight Wij . The connectivity of the network is repre- 
sented by the adjacency matrix A whose element Aij = 1 
if there is a link from a node j to i. Otherwise, Aij = 0. 
We set An = conventionally. The degree ki of a node i 
is given as ki = ■ Aji . Since we consider weighted net- 
works with weight Wq. we define the weighted adjacency 
matrix A as Aij = Aij . 

The motion of a walker on the weighted networks de- 
fined by the matrix A is a stochastic process in the dis- 
crete time. We derive the stationary distribution P°° of a 
walker being at node i following the method of Ref. [24| • 
To set up the equation, we define the transition probabil- 
ity as follows. A walker at node i at time t selects one of 
its ki linked nodes with hopping probability X^. Then, 
at time t + 1, the walker moves to the selected node. The 
hopping probability Tji from node i to j is then given as 
Tji = Aji j Ki . where Ki = J^j Aji is the strength of node 
i. As an initial condition, assume that the walker starts 
at the node q at time t = 0. Then the recurrence relation 
of the transition probability P, g of finding the walker at 
node i at time t is 

P iq (t + l) = J2 T U P k( t ) ■ ( 2 ) 
l 

Then the transition probability Pi q (t) is written by iter- 
ating as 

P iq (t)= E T a t - l ---T hh T hq . (3) 

For a symmetric A with A^ = Aji, one finds K q Pi q (t) = 
KiP q i{t) by comparing P q i and Pi q . In the stationary 
state, the probability P°° of finding a walker at node i 
should be independent of initial starting nodes, which 
gives KiP^° = K q P?°. Summing up over q, one finds 

pr = k^n , (4) 

where N = J2 q =i K q = Y. q =i ELi In wei g hted 

networks with symmetric weights, P°° is proportional to 
the strength of node i, i.e. the sum of the weights of the 
nearest neighboring nodes. The same result was found 
in the recent study on the dynamics of random walks on 
growing weighted networks [22| . 

In this paper, we consider the symmetric weight Wy, 

Wij = (kikj) a . (5) 
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For the weight P°° is not given as a simple form. 
Hence it is better to handle the distribution P^° of finding 
a walker at nodes with degree k. Using Eq. (j4]), one can 
see that 



poo 



JV 

£ 

i=l 



1 



N N 



^^A^hk.rS^k (6) 

i=i i=i 



To express the sum in Eq. ([6]) in terms of degree fc, we 
arrange the sum as follows. Only terms with ki = k 
contributes nontrivially to the sum and thus NP(k) 
nodes with the degree k in a network have the nontriv- 
ial contributions to the sum. The node with the degree 
k has k linked neighbors whose degree ranges from 1 to 
the maximal degree of the network k max . Hence, the 
number of nontrivial terms in the sum Yli=i ^2j=i * s 
NP(k)k, which can be arranged in the order of increas- 
ing degree. Then, the double sum of Eq. is written as 
NP(k)k a+1 (. g (l)l +g(2)2 a + ■ ■ • where 
g(k') is the degree distribution of the node involved in 
such NP(k)k terms. For large N, we approximate g(k') 
to P(k'). Then P£° is approximately given as 



PS 



JjP{k)k a + 1 J P{k')k' a dk' 



P(k)k a+L / 



P{k')k' a+1 dk' . (7) 



On SFNs with degree distribution P(k) ~ k 7 , the inte- 
gral in the second line is finite for a < 7 — 2. Hence we 
finally obtain P fc °° on WSFNs as 



Pi 



k" 



1 



7 



(8) 



The exponent a varies with a and 7, and also changes its 
sign. For a = 7— 1, i.e. a = 0, Pj^° is independent of de- 
gree k so a walker does not feel the inhomogeneity of the 
underlying network structure. While a walker performs 
biased walks to nodes with the larger degree for a > 0, 
the direction of the bias is reversed for a < 0. Since 
the exponent a is a free parameter, one can controls the 
direction of the bias for a given 7. 

To check the scaling relation |(5J), we perform Monte 
Carlo simulations on the WSFNs with N = 10 5 and the 
average degree < k >= 4. In the steady states, we mea- 
sure P£° for various a up to 2.6 for 7 = 2.7 and 3.0 for 
7 = 3.3. Fig. 3 shows the plot of against k for sev- 
eral values of a. As shown, P£° scales in power-law with 
k. The inset in each panel shows the plot of a against 
a. The simulation results agree well with the analytical 
prediction ©. 
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FIG. 3: (Color online) The plot of P fe °° and a for 7 = 2.7 (a) 
and 7 = 3.3 (b). Insets show the relation |(8} (solid line) and 
numerical estimates of a (symbols). 



weights ([5]). A lamb and a lion initially locate separately 
on randomly selected two nodes. Then the probability 
D(k) of finding two walkers at the same node with de- 
gree k at the same time is proportional to (P£°) 2 . From 
Eq. ©, one gets 



with 



D(k) = (pn 2 /NP(k) ~ k v 



v = 2(a+ 1) -7 



(9) 



(10) 



Then the probability D of finding two walkers on the 
same node with any degree is given as 



D = 



D{k)dk 



k 2(a+l)- ldk (11) 



fco 



Since the upper bound k max diverges with TV, the integral 
jkma* k 2(a+i)- ldk diverges for a > (7 - 3)/2. Hence 
there exists a crossover value a c given as 



IV. CAPTURE PROCESS ON WSFNS WITH 
SYMMETRIC WEIGHTS 

In this section, wc consider the capture process or 
the lamb-lion problem on WSFNs with the symmetric 



(7 - 3)/2 



(12) 



For a < a c , the lamb survives indefinitely with a finite 
probability. However, for a > a c , the lion captures the 
lamb with the unit probability. To check the scaling rela- 
tion fTU]). we measure D(k) on the WSFNs with N = 10 5 
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FIG. 4: (Color online) The plot of D(k) and v for 7 = 2.7 (a) 
and 7 = 3.3 (b). Insets show the relation of (|10p (solid line) 
and numerical estimates of v (symbols). 

and < k >= 4. For 10 5 trials, we count the number n(k) 
of capture events on nodes with the degree k. We obtain 
D(k) to divide n(k) by total trials(10 5 ). Fig. 4 shows 
the plot of D(k) for several values of a, which scales well 
with k in power-law. As shown in the insets of Fig. 4, 
numerical estimates for v satisfy the relation (|10[) very 
well. 

To verify the existence of a c by another method, we 
now consider the survival probability S(t) of a lamb. S(t) 
always satisfies S(t) — e~*/ T on random and scale- 
free networks due to the small world nature [H, HH ■ As 
S(t) = Soo e~ t / T in SFNs with any 7, we are interested in 
the average life time (T) of a lamb rather than S(t) itself. 
From (T) = f™[-dS(t)/dt] dt and S(t) = 5 m e~*/ r , we 
have (T) ~ r. Hence (T) is infinite for a < a c and finite 
for a > a c in the limit N — > 00. However, for the finite- 
sized networks, a lamb is eventually captured within N 
time steps for any a. For a < a c , the maximum life time 
should be the order of N to guarantee the finite survival 
probability in the limit N — * 00. Hence (T) is expected 
to scale as (T) ~ N for a < a c . We measure (T) on 
WSFNs of 7 = 2.7 and 3.3 with network size N up to 
10 6 . From Eq. (T5J), one reads q c = -0.15 for 7 = 2.7 
and a c = 0.15 for 7 = 3.3. 

In Fig. 5, we plot (T) against N. As shown in each in- 
set, S(t) exponentially decays for any a. For 7 = 3.3(Fig. 



FIG. 5: (Color online) The average lifetime (T) of a lamb for 
7 = 3.3 (a) and 7 = 2.7 (b). The solid line is a guide to 
the eye. Insets show the semi-logarithmic plots of S(t) for 
TV = 10 s . From top to bottom, each line corresponds to the 
S(t) of a < a c , a c and a > a c respectively. 



5(a)), (T) increases with N as with <j) = 0.94(1) for 
a = 0.05 (< a c ) and = 0.90(1) for a = a c {= 0.15). We 
estimate <j> by measuring successive slopes from the log- 
log data in Fig. 5(a). For a = 1.0 (> a c ), (T) tends to 
saturate to the asymptotic value (T)oo with decreasing 
successive slopes. The exponent <j> of a = 0.05 is close 
to the expected value 0=1. For a = a c (= 0.15), (T) 
seems to diverge with <fi = 0.9. However, since (T) of 
a = 1.0 already tends to saturate, it is expected that (T) 
of a > a c would saturate to a finite value in the network 
with N > N c (a), where N c (a) is the characteristic size 
for given a. For example, for a — 1.0, (T) does not get 
into the saturation region even after N = 10 6 , which im- 
plies -ZV C (1.0) > 10 6 for a = 1.0. Since N c should increase 
as a — * a c , it is empirically impossible to see the satura- 
tion of (T) via simulations. Therefore, the initial slope 
4> at a c may have no special meaning as that of a > a c . 
The same behavior for (T) was observed for a = case 
[lH , where (T) initially algebraically increases with con- 
tinuously varying (/>(< 1) as 7—5-3 from below. 

For 7 = 2.7 (Fig. 5(b)), we estimate 4> = 1.00(2) for 
q = —1.0 (< a c ) as expected. However, for a = —0.05 (> 
a c ), (T) algebraically increases with <\> = 0.82(1). Since 
for a = N c is already larger than 10 6 for 7 = 2.75, 
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• a = 0.05 

a= 1.0 




FIG. 6: The plot of m(k) for 7 = 3.3. The solid and the 
dashed line correspond to a = 1.0 and 0.05 respectively. The 
inset shows the scaling plot m(k)k~ v with r\ — 1.06 for a 
0.0") (dashed line) and r) — 1.95 for a = 1.0 (solid line). 



it is difficult to see the saturation of (T). For a = a c , we 
estimate cj) = 0.90(1). As in 7 = 3.3, the initial slope for 
a > a c has no special meaning. Based on our numerical 
results, we are convinced that (T) approaches a finite 
value (T^oo for a > a c and becomes infinite for a < a c 
in the limit N — ► 00. Hence in the limit N — > 00, we 
have 



lim S(N,t) 

N->oo 



' S e- f / T - (a>a c ) 
(a < a c ) 



S Q 



(13) 



with the characteristic time 7% 



V. AVERAGE MASS OF A NODE WITH 
DEGREE k 

Another interesting quantity in condensation phenom- 
ena on networks is the average mass m(k) at a node with 
degree k in the steady state [H, [H, Ell ■ In ZRP with 
chipping rate u(m) ~ m s , the complete condensation 
takes place for S < S c , where S c = 1/(7—1) for un- 
weighted SFNs fijl and' (a + l)/( 7 - 1) for WSFNs with 
the weight (0 pj|. For S < S c , m(k) increases as k a+1 
for k < fc c , and as k^ a+1 ^ s for k > k c on the WSFNs. 
Especially for S — 0, m(fc) increases k a+1 until /c < k max 
and jumps to the value rrihub ~ p^V at k max - Hence the 
condensation takes place at the node with k max degree 
in ZRP. 

The recent study on CA model on unweighted SFNs 
showed that m(k) linearly increases up to k ma x without 
the jump at k max unlike in ZRP with constant chipping 
rate The linearity of m(k) comes from the fact that 
all masses can diffuse. The mass rrihub formed on the 



node with the degree k max can diffuse throughout net- 
work to make the steady sate distribution P?° . By taking 
average over all nodes, the rrihub soaks into the average 
mass in(k) unlike in ZRP where all samples have irihub 

at 

kmax- The linearity of m(k) on unweighted SFNs results 
from P~ - kP{k) pi]. Therefore, from P fc °° ~ k a+1 P(k) 
on the WSFNs, we expect m(k) ~ k a+1 up to khub- To 
see this explicitly, we derive the relation ~ k a+1 as 
follows. 

We consider the average total mass M(k) of nodes with 
degree k defined as 



M(k) = mP^m.k) , 



(14) 



m=0 



where P 00 (m,k) is the probability of finding a walker 
with mass m at nodes with degree k in the steady state. 
Since the mass distribution P(m) in the steady state is 
independent of k, we write Poo(m, fc) = P(m)P£°. From 
([7]) and (T14"]) , one reads 



M(k) 



k a+1 P{k) Y mP(m) , 

m=0 



(15) 



where we drop the normalization constant of P^° . Since 
the number of nodes with degree k is NP(k), m(k) is 
given as 



m(k) 



M(k) 
NP(k) 



(16) 



To confirm the scaling behavior of m(k), we measure 
m(k) in the condensed phase on the WSFNs of 7 = 3.3 
and N = 10 5 . In Fig. 6, we plot m(k) against k for 
a = 0.05 and 1.0. With w = 1, we set p = 3.0 which 
corresponds to the condensed phase for both a values. 
Assuming m(k) ~ k v , one expects r\ = 1.05 for a = 0.05 
and 2.0 for a = 1.0 respectively. From the scaling plot 
m(k)/k T ' (Inset of Fig. 6), we estimate rj = 1.06(2) for 
a = 0.05 and 77 = 1.95(5) for a = 1.0 which agree well 
with the predictions. 



VI. SUMMARY 

In summary, we investigate the properties of 
conserved-mass aggregation (CA) model on weighted 
scale- free networks (WSFNs). In WSFNs, the weight 
Wij is assigned to the link between node i and j. We 
consider the symmetric weight given as Wij = (kikj) a . 
In CA model, masses diffuse with unit rate and unit 
mass chips off from mass with rate ui. In addition, the 
hopping probability Tji from node i to j is given as 
T oi =Wji/T, <m> w mi . 

On the WSFNs, a walker finally reaches the hub node 
with k max degree for a — > 00, while it is trapped forever 
at nodes with the minimal degree for a — > — 00. In the 
lamb- lion capture process, it means that the lion captures 
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the lamb at the hub node within finite time interval for 
a — ► oo. On the other hand, a lamb survives indefinitely 
with finite probability for a — > — oo. because the lion 
cannot escape from a node with the minimal degree to 
capture a lamb at some other node. In-between the two 
limits, one expects a crossover a c below which the life 
time (T) of a lamb is infinite. However, for a > a c , (T) 
is finite. The dependence of (T) on a is similar to that 
on unweighted SFNs of a = where (T) is infinite for 
7 > 3 and finite otherwise [l3[ . 

To verify the existence of a c , we need the stationary 
distribution of finding a walker at nodes with degree 
k. From the equation for the transition probability Pji(t) 
to go from node i to j in t time steps, we analytically find 
P£° ~ fc Q+1 ~ 7 . Next, we consider the so-called lamb-lion 
capture process. With we find the probability D{k) 
of finding two walkers at the same node with degree k 
at the same time to scale D(k) ~ fc 2 ( Q+1 )~ 7 . Finally, 
integrating out D(k), we find the death probability D 
of a lamb. A lamb survives indefinitely with the finite 
survival probability for a < a c , while it is eventually 
captured by a lion for a > a c . We analytically find 
a c = (7 — 3)/2. Therefore, in the limit N — > 00, the life 
time (T) of a lamb is finite for a > a c , while it is infinite 
for a < a c . We numerically confirm the all analytical 
results. 

The existence of the condensation transitions is known 
to depend on (T) of a lamb For a > a c , (T) is 

finite so the condensation always occurs for any nonzero 
density. On the other hand, for a < a c , the infinite (T) 



ensures the condensation transitions at a certain critical 
density p c . For a > a c , we numerically confirm that the 
condensation always takes place at very low density. We 
also numerically confirm that for a < a c , CA model on 
the WSFNs undergoes the same type of the condensation 
transitions as those of SCA model in regular lattice. 

Finally, we investigate the behavior of the average mass 
m(k) of a node with degree k. In ZRP with constant 
chopping rate on networks [HI, [l6j|, m(k) increases as 
fc Q+1 , and jumps to the total mass of the system at k max . 
However, in the SCA model on unweighted SFNs, it was 
shown that m(k) linearly increases with k up to k max 
without any jumps fl3l | . Furthermore the linearity of 
m(k) is valid for any p > 0, which comes from the fact 
that the diffusion is only the relevant physical factor to 
decide the distribution m(k). Similarly, on the WSFNs, 
we analytically find and numerically confirm that m(k) 
algebraically increases as k a+1 for any p > without any 
jumps. 
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